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ABSTRACT 



D ' Aims. In order to understand the anisotropic properties of local radiation field in the curved spacetime around a rotating black hole, 

' we investigate the appearance of a black hole seen by an observer located near the black hole. When the black hole is in front of a 

source of illumination the black hole cast shadow in the illumination. Accordingly, the appearance of the black hole is called the black 
hole shadow. 

Methods. We first analytically describe the shape of the shadow in terms of constants of motion for a photon seen by the observer 
in the locally non-rotating reference frame (LNRF). Then, we newly derive the useful equation for the solid angle of the shadow. In a 
[•j'] [ third step, we can easily plot the apparent image of the black hole shadow. Finally, we also calculate the ratio of the photon trapped 

by the hole and the escape photon to the distant region for photons emitted near the black hole. 

Results. From the shape and the size of the black hole shadow, we can understand the signatures of the curved spacetime; i.e., the 
(-H ' mass and spin of the black hole. Our equations for the solid angle of the shadow has technical advantages in calculating the photon 

' trapping ratio. That is, this equation is computationally very easy, and gives extremely precise results. This is because this equation 

I , is described by the one-parameter integration with given values of the spin and location for the black hole considered. After this, the 

Q ■ solid angle can be obtained without numerical calculations of the null geodesies for photons. 
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1. Introduction 

K> To understand the general relativistic effects on the radiation field around a black hole it is essentially important to elucidate the 
, nature of both the physics and the astrophysics in the curved spacetime. This is because the radiation field plays energetically 
■ and dynamically important roles in the accretion flows plunging into the black hole, especially for the supercritical or hypercrit- 
ical accretion flows, where the photons or neutrinos interact with matter. Even for black holes with low mass accretion rates, the 
general relativistic effects of the radiation field becomes important because the observational signatures directly reflect the nature 
of the radiation field. So far, the radiative transfer in curved spacetime was investigated by many authors (e.g. Lindquist 19661 
O . Anderson& Spiegel [19721 S chmi d-Burg k [19781 T horne [19811 S chinder [19881 Tu rolla & No bili.1988; Anile & Romano , 19921 
O . Cardall & Mezzacappa [2003l Park [2006l Takahashi [2007bl [20081 De Vilfiers [20081 Farris et al. [2008l l. 

I ' It is believed that the spacetime geometry around a black hole in the real universe is well described by the Kerr metric. In this 
^ , case, clarifying the effects of the black hole spin on the radiation field is required to understand the nature of the radiation field in 
the black hole spacetime. Any radiation field around the rotating black hole is influenced by frame-dragging effects, which make an 
anisotropy of the local radiation field. That is, the anisotropic features of the local radiation field directly reflects the effects of the 
^ black hole spin. In the optically thick limit of the radiation field, it is expected that the anisotropy becomes very small (but not zero), 
■ ■ because the radiation field and the matter intensely interact. On the other hand, in the optically thin limit of the radiation field the 
anisotropic nature can be clearly seen in the vicinity of the black hole. For example, for the black holes in the center of the gamma- 
ray bursts (GRBs) the effects of the neutrino annihilations, which are potentially the energy source of the explosion of the GRB, are 
well investigated in the optically thin limit of the neutrino radiation field (e.g. Asano & Fukuyama'2001', Miller et al. '2003' Birk et 
al.|2007). In this case, the black hole's rotation enhances the deposition energy due to the neutrino annihilations. In the optically thin 
limit, the photons or neutrinos emitted at some location near the black hole are transfered along the null geodesies and some of them 
are swallowed by the black hole. As is well known, because the propagation rays are influenced by the frame-dragging effects due 
to the black hole's rotation, at some location near the black hole the effects of the black hole spins are directly seen in the shape and 
the size of the appearance of the black hole. When a black hole is in front of a source of illumination, the black hole casts a shadow 
on the illumination. Accordingly, the appearance of the black hole is sometimes called a black hole shadow (e.g. Falcke, Melia & 
Agol 2000; Takahashi 2004, 2005; Huang et al. 2007; Bambi & Freese 2009; Hioki & Maeda |2009l or silhouette (e.g. Fukue l2003l 
Broderick & Loeb 2009; Schee & Stuchlik ,2009aJ . With regard to the realistic accretion flows around black holes, when we observe 
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the black hole with accreting gases, we only see the escape photons emitted by surrounding gases, while the trapped photons make 
up the shade region. 

The observers near the black hole see the black hole shadow in their sky. In this paper we calculate the analytic expression of 
the shape of the shadow seen by the observer located near the black hole. We also obtain the new equation for the solid angle of the 
shadow in the sky of the observer. The main purpose of the present study is to show the analytic expressions of the shadow and the 
equation for the solid angle. Our equation for the solid angle of the shadow has technical advantages. By using this equation, the 
calculations of the null geodesies around the Kerr spacetime is not required for the calculations of the solid angle. Also, we do not 
need the particles to be emitted and judge which particles will be swallowed by the black hole, as in the calculations done in e.g. 
Thorne (11974 ). Since the equation for the solid angle presented in this paper is described by the one-parameter integration for the 
given values of the black hole spin and the location, this equation is computationally very easy. We assume c = 1 and G = 1 . After 
giving some preliminary calculations in Sect. |2] in Sect. |3]we give the analytic expressions for the shape of the shadow seen by the 
observer near the black hole. In Sect.|4] the equation for the solid angle of the shadow is presented. We give the conclusions in Sect. 

s 



2. Preliminaries 

2.1. Background metric 

We describe the background metric for a Kerr black hole by using the spherical Boyer-Lindquist coordinate. The 3+1 form of the 
metric g^y is described as 

ds^ = -a^dP- + yijidx' + 0dt){dx^ + /3^dt). (1) 
On the other hand, the inverse of the metric g''''' is given as 

g'^^'d.dy = ^(d,- 0di){d, - p^dj) + yjdidj. (2) 

For the Kerr metric described by the spherical Boyer-Lindquist coordinate, the non-zero components of the lapse function (a), the 
shift vector (j3' and y6,) and the spatial components of the metric and y'-') are given as 



A A p^ p^ Aiiivte 

and jij - yj - for all other /, j. Here, A - - 2Mr + a^, p^ - + cos^ 6, A - (r^ + a^)^ - a^A sin^ and oj - 2Mar/A and 
we use yS; = jij/S-' and y'^ytj = 6'j. The determinant of y,j is given as dety,y = yn-yeey^^ - (y'''''y'"'y'*'*)"' = Ap^ sin^ 6/ A. 



2.2. Locally non-rotating reference frame tetrads 

The locally non-rotating reference frame (LNRF) is the frame moving with the four velocity given as - ad'f. For this frame, the 
tetrad vectors are given as 

ef^ = [a, 0, 0, 0] , = [0, Vtv, 0, 0] , = [0, 0, Vt^, 0] , = Vr^[y8*, 0, 0, l] . (4) 

In the same way, we have 

e'^o)--[h 0, 0, -yS^] , e^'i) = [o, Vr^, 0, o] , 4, = [o, 0, V7^, o] , 4) = [O, 0, 0, Vr*] . (5) 



2.3. Photon four momentum 

The photon four momentum of the photon is calculated from the action S as pa - daS . The action S is given as 



1 vR r 

S = -m^A-Et + L-if> + s,- j —dr+sg\ ^/®de. 



(6) 



where A is the afRne parameter, m is the mass of the particle, E is the energy and is the angular momentum with respect to z-axis. 
The signs Sr{- ±1) and sg(= ±1) correspond to the directions of the propagation of the null geodesies, i.e. the signs of pr and pe 
determine the signs Sr and sq. The functions &(ff) and R{r) are given as 

0(0) = Q-cos^e[-a^E^ + Ll/sm^6), (7) 

R(r) = [(r^ + a^)E - aL,f - A [{L, - aEf + Q\ . (8) 

Here, Q is the Carter constant. Along the null geodesies, E, L, and Q are constants of motion. The photon momentum is explicitly 
given as 

I- 

p, - d,S - -E, p^ - d^S - Lj, Pr - drS - ^r—, Pe - deS - sg V©. (9) 
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Hereafter we newly define 



n = —, a, = —, A. = — = r: - Ir, + a' A* = —, to, = ojM, pi = (10) 
and 

then we obtain 

R.{r,) = [{rl + flj) - - A. [(^ - a^^f + //] , (12) 
0,(0) = T] + (a,-O^-{asm0-(/smef. (13) 

Since cannot be negative, there is a constraint of ?; + (a, - ^)^ > 0. Now, for 77 - (a, - ^)^ - 0, we have sin^ — (.la^; that is, 
6 =constant. Hereafter, we drop the asterisk of r» for simplicity. The photon four momentum in the LNRF is calculated as 



2.4. Innermost unstable orbit 

In the motion of the r-direction, the radius of the innermost unstable orbit rj for a photon is calculated from the equations given as 
Rti^s) = and drR»{rs) - 0. From these equations, we have 

^^'IzAzll^^^ rj^ [4AJr,)-rs(r,-lf] for a. ^ 0, 7, + ^ 27, r = 3 for a, = 0. (15) 

These results are presented in the past studies (e.g., Bardeen|T973j Chandrasekhar |T983l l . We note that the innermost unstable orbit 
was calculated in a different way in the recent study by Schee & Stuchlik ( |2009al l,Fl in which deeper explanations of the innermost 
unstable orbit are given (see Sect. 3.3 in Schee & Stuchfk l2009al l. 

3. Shape of the black hole shadow as seen by the locally non-rotating reference frame observer 

Here we calculate the shape of the black hole shadow in the sky of the LNRF observer near the black hole. In order to do this, we 
first introduce the angular coordinates (0, 0) in the sky of the observer. The direction of the null geodesies passing the observer is 
described by these coordinates. That is, the photon momentum in the LNRF can be expressed as 

p'^^^E, = cos 0, = £■ sing cos 0, = £ sin sin 0. (17) 

These relations define the angular coordinates (6, (p). In these definitions, the direction of 5 = corresponds to the direction of the 
black hole. Here, we assume 



,.(0) ^(1) ^(2) ^(3)-, 



<p]. (18) 



In this case, the direction of p^''^ corresponds to the angle 6 - tt and the direction of p*^* corresponds to = 0. 
On the other hand, by using the LNRF tetrads given above, we have 



A, sin0^ 



1 - w*^ 



The parameters ((, rf) describing the shape of the black hole shadow are given by Eq. ( fTSl l (see, e.g. Bardeen[T973j Chandrasekhar 
I1983I I. From the definitions of {9, (p) in Eq. (fTTl l. we have 



g = COS-' , ^ = tan-' ^ 



pie) /pit) 



(20) 



' In Schee & Stuchlik (2009a'l, the parameter £,(= + rf) is used as one of the constants along the geodesies. By using this, the functions R,(r,) 
and dr.Rtir,) are related to X and d^.-C as 

£ - 2a,( -a, + , dr,£ = d,-,R, + [2r,A, - (r, - l)(r; + a, - fl.^)J . (16) 

By using R,(r,) = 0, d,-,R(r,) = and ( given in Eq. dlSb . it is shown that -Cmax = and c?,-.-£max = 0. Here, when R,(r,) = 0, .£ = -£max where 
Xmax is given by Eq. (41) in Schee & Stuchlik (i2009aj . From these calculations we can see that the calculation method for the innermost unstable 
orbits used in Schee & Stuchlik ( i2009ab produces the same results in terms of ^ and 77 as in the past studies of Bardeen ( 11973b and Chandrasekhar 
([1983}. 
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Fig. 1. Shapes of the black hole shadow seen by the observer near the black hole. In this plot, the normalized Lambert map 
projection described by Eq (l24l i is used. The abscissa is 9^ and the ordinate is 0,.. In this figure, the direction of the black hole (9 - n) 
corresponds to (0^, 9y) - (0, 0), and the opposite direction (9 - 0) corresponds to the outermost circle (i.e., 0^ + 0^ = 1) denoted by 
the dotted line. The projection of the rotation axis of a black hole corresponds to 9jc = 0. The locations of the observers are 9 = 85° 
and r - 100, 10, 3, 2, 1.3, 1.1, 1.05, 1.045. The black hole spin is a* - 0.999. In this case, the radius of the event horizon is 
~ 1.0447M. 



With these relations we can specify the position in the sky of the observer from the photon momentum and vice versa. By using Eqs. 
( fT9] ) and ( l20l i. we can draw the contour of the black hole shadow in the sky of the observer located at (r, 9). Along the contour of the 
shadow, we change the parameter of Eq. ( fTsT i in the range of ©» > (see, e.g. Bardeen [19731 Chandrasekhar [1983]) . Hereafter, 
we introduce new coordinates 

X - sin9sm{-<p), F = sin 0cos(-^), Z-cos9. (21) 

With the coordinates (X, Y,Z), we plot the projection of the black hole shadow onto the sky. This contour is expressed by a closed 
curve in the sphere (see Fig|2j discussed below). Before doing this, we have to specify the signs of Sr and sg. To plot the contour of 
the shadow, there are two angles of @(9) - 0, where the sign of sg changes across these angles. When the observer is located near 
the black hole, there are points where the sign of Sr changes. The sign of sg is determined by 

ifl = +1 for Y>0, for F < 0. (22) 

On the other hand, when the observer is located far from the black hole, we have Sr = - 1 . This means that all the photons swallowed 
by the black hole move toward the black hole, i.e. p'' < 0. However, when the observer is located near the black hole, even if the 
photons are emitted in the outgoing direction p'' > 0, some of them are swallowed by the black hole. In this case, there are points 
where the sign of Sr changes in the contours of the shadow. The sign of s,- is determined as 

= +1 for drR{r) > 0, = -1 for dr R(r) < 0. (23) 

We note that in Rir), only ( and rj depend on the value of r,. 

The three-dimentional description {X, Y, Z) can be projected with the Lambert map projection (see, e.g. Oihe l2000] l given as 



0v = fccosMsinv , , / 2 , ( m = sin F 

n I ■ here k= a/- , and < -i,,, ,™ ■ (24) 

9y^ksmu Vl+cosMCOsv [v^tan^X/Z) 

In Fig.dl we plot the shapes of the shadows with the normalized Lambert map projection. The abscissa is 0i and the ordinate is 9y. 
In this projection given by Eq (l24l i. all the area of the sky of the observer is mapped into the region of 9^ + 9^ < 1 . In this figure, the 
direction of the black hole (9 = tt) corresponds to (9^, 9y) = (0, 0), and the opposite direction (9 = 0) corresponds to the outermost 
circle (i.e., 9^^ + 9^. - 1) denoted by the dotted line. Both 9^ and 9y are normalized so that their maximum values become unity. The 
locationoftheobserveris6' = 85°andr = 100, 10, 3, 2, 1.3, 1.1, 1.05, 1.045. The black hole spin is a. = 0.999. In this case, the 
radius of the event horizon is r+ ~ 1 .0447M. In this figure, the outermost circle denoted by the dotted line is the line of 9 - 0. As 



5 






Fig. 2. Black hole shadow in the sky of the observer near the black hole in the three dimensional plot (X, Y, Z). The parameters 
are the same as those in Fig.[T] The position of the observer is at the origin {X, Y,Z) = (0,0,0) (filled triangles). In the sky of the 
observer, the direction of the black hole is (X, Y,Z) - (0,0, 1). The shaded regions are shadows for r = 10. The contours of the 
shadows in the front (opposite) side of the sphere are plotted by the sohd (dashed) lines. 
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expected, the size of the shadow becomes larger for the observer located closer to the black hole. Interestingly, we can see that some 
photons in the region of 0^ < 0, can escape to infinity even when the observer is located very near the black hole. This is because of 
the frame-dragging effect. It is confirmed that no photons emitted at the horizon can escape. In Fig.|2]we plot the projection of the 
shadows onto the sky in the three dimensional coordinate {X, Y, Z), where the observer is located at the origin, (X, Y, Z) - (0, 0, 0) 
in this plot (filled triangles), and the direction of the black hole is {X, Y, Z) - (0, 0, 1). While most of the photons emitted from the 
location far from the hole (e.g. r = 100) escape to infinity, many of the photons emitted near the black hole are swallowed by the 
hole (e.g. r - 1.05 and 1.045). The parameters are the same as in Fig.[T] The similar signatures seen in Fig.[T];an also be seen in this 
figure. In Fig|2] the shadow has a solid angle and its shape is deformed like in Fig.[T] In the next section, we give the new equation 
for the solid angle. 

4. Solid angle for the photons trapped by the black hole 

Here we give the equation for the solid angle of the black hole shadow in the sky of the LNRF observer. This solid angle is the 
function of the black hole spin and the location of the observer (r, 0). When the observer is moving with some velocity, the solid 
angle also depends on the velocity. Here we only consider the LNRF observer for simplicity. Below we derive first the new equation 
for the solid angle of the black hole shadow in the sky seen from an observer located at a point near the black hole. Next, we show 
the distribution of the ratio of the solid angle in the sky in the x-z plane. Note that this ratio corresponds to the ratio of the trapped 
photon by the hole and the escape photon toward distant regions, when we consider the photons isotropically emitted from a point 
(as a radiation source in LNRF) near the black hole. 



4.1. Equation for the solid angle 

The shape of the black hole shadow is symmetric with respect to the Y-Z plane. This makes it sufficient to calculate the solid angles 
surrounded by the 0i-axis and the contours of the black hole in the region of F > 0. The twice the amount of the solid angle for 
7 > corresponds to the amount of the solid angle for the whole region of the black hole shadow. First, we consider the case when 
the whole region of the black hole shadow exists in the region of Z > 0, which is the front side of the observer when the observer 
faces the black hole. In this case, we have if = -1 and the equations describing the sphere in Z > are given as 

Z^(l-X'-Y'y''.r(X,Y), f^^^^—±-^^, /;.__ZI__, (25) 

where fx and fy are the partial derivative of the function / with respect to X and Y, respectively. Now, the solid angle S of the 
region Q of the shadow on this sphere is calculated as 

Here we introduce the variables (r, 6) as 

X^rcosO, Y^rsinO. (27) 
From these definitions, we have 



r-5 r - Y -1 de dd> 

r = Vz2 + y2 ^ sine, tan0=- = ^, — = 

X tan(/> ors ors 

By using these, the solid angle S is calculated as 

S ^2 f f ^2 f f , drdf) = 2 f d0 f ' , ^ dr = 2 f dell - Jl^, 
J Jn J Jn Vl Jo Jo Vl -r2 Jo L 

Because we now consider the case of Z > 0, i.e. < < 7r/2, we obtain 

-y/i"^^^^^^^ = Vl - sin^e = Vcos^e = cos §. 
Next, by changing the integration variable from to r,, we have the ranges for the integration as follows 

61 = ^ TT and = r™" ^ rf", with t/fl = t/r,. (31) 
Then, the solid angle S is given as 

5=2 Jj^ dr, (-|^ W - cos 0), (32) 



(28) 



(29) 



(30) 



where 

/ 1 \'^' 

COS0 ^-sA — \ (here = -1). (33) 

\A,) 1 - w,^ 
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Now we we consider the case where the whole of the black hole shadow exists in the region of Z < (the back side of the 
observer). In this case, we have s,- = - I and the equations describing the sphere in Z < are given as 

Z^-il-X'-Y')'''^f-iX,Y), f^^-—l—^^, = (34) 
In the same way, we have 



- '"max(S) = Vl -Sin^e = VcOs2 = - COS 0. (35) 

The change of the sign of do not cause any change because both and fy are used in the form of (/^)^ and {fy)^, respectively. 
In terms of the change of the sign before cos 9, this change is canceled out by the change of the sign of Sr- Consequently for the 
shadow in Z < the equation is essentially same as for the shadow in Z > 0. For the shadow crossing the Z - Q, the equation 
assumes the same form. 

In summary, the equation of the solid angle S is given by 



c/rj (-5^,0) (1 -cos 0), where cos = -i,. I , ^* . iwcj)^—^^ — (— P= I • (36) 

\A,/ Ay sinflWC/ 



Here the range of the integration for is determined by the condition of 0, > 0, which can be numerically calculated. Moreover, 
by using the analytical form for ^, rj, 0,, ^, and so on, the integrand of the integral of Eq. ( |36] ) is analytically given as 



S = -4p^sin6i J^^ 



irl - 3r2 + 3r, - al)[rl{r, - 3)^ - fl,(r, - \){r^ + al cos^ 0)] (^^^^ ^ 



ai(rs - l)HA sin^ 00. + p4^2) a/©! \ 1 - w.^ 

The derivation of Eq. (l37b is given in App.lAl The Eqs. (|36] | and (|37] | are the main results of the present study. 



dr.,. (37) 



4.2. Photon trapped ratio 

Based on Eqs. ( |36] | and (l37T l. we numerically calculate the distribution of the photon trapped ratio S /{4n) around the black hole. In 
Figl3]we plot the trapped ratios S /{4n) for a, = (top) and 0.999 (bottom) in the x-z plane. For the rapidly rotating black hole, the 
photon trapping ratio decreases especially near the equatorial plane because of the frame-dragging effects. That is, the photons near 
the equatorial plane can escape to infinity more easily than the photons near the pole. These features can also be seen in the shapes 
of the shadows. Surprisingly, as shown in these plots, 10 % to 30 % of the photons emitted very close to the black hole (e.g. r < 1.2 
for fl» = 0.999) can escape to infinity. Although these calculations do not take into account the effects of the gravitational redshift, 
these escaping photons can transport the information of the spacetime and the general relativistic phenomena to the observer at 
infinity. These results about the trapped ratio for photons are very important when we consider the radiative phenomena that arise 
near the event horizon. In Fig.|4]we plot the lapse function a representing the effects of gravitational redshift for black hole spins 
of a/M = (top) and 0.9 (bottom). The contours of the lapse function in Fig. |4] shows nearly spherical shapes, which is contrast to 
the plots of the photon trapped ratios in Fig. [3] Then, even if we take the effects of the gravitational redshift into account, the energy 
amount of the escaping photons near the equatorial plane is than that near the pole region for the case of the rapidly rotating black 
hole. 



5. Discussion and conclusions 

Since the last century, observational data containing information of the black hole spacetime have been obtained for the black 
hole candidates in e.g. the active galactic nuclei, the X-ray black-hole binaries and the Galactic center of our galaxy. However, the 
parameters of the black hole spacetime like mass and spin estimated in the past studies usually and highly depend on the assumed 
accretion flow models, and so far the final values of the parameters of the black hole, especially the black hole spin, have not been 
obtained for any black hole candidates. This is because the observational data usually contain the information of the plasma physics 
of the accreting matter in addition to the information of the spacetime metric, and frequently both these informations (metric and 
plasma physics) are degenerate in the observational data. In these cases we cannot uniquely determine the spacetime metric from 
the observational data. Because the optical phenomena directly reflect the spacetime metric, the spacetime metric can be uniquely 
determined if the observational signatures closely related to the optical phenomena like the contours of the black hole shadow can be 
found in the observational data. For this purpose a deep understanding of the optical phenomena around the black hole is important. 

In past studies it was proposed that the black hole's metric can be estimated from the observational features like the line and 
continuum spectrum, quasi-periodic oscillation (QPO), radio visibility, polarization (see, reviews, e.g., Remillard & McClintock 
120061 Psaltis 2008). However, all these observational signatures highly depend on the accretion flow models. Actually, for some 
black hole candidates, several values of the black hole spin can be proposed for one and the same black hole candidate. For example, 
in terms of the iron line spectrum observed in an X-ray black hole binary, it is pointed out that the Fe K line profile and the estimated 
value of the black hole spin of the black hole binary GX 339-4 highly depend on the assumed continuum spectrum, which is first 
subtracted from the observational data to make the line profile (Yamada et al. 2009). This means that the black hole spin cannot be 
correctly estimated from the line profile unless we can uniquely determine the accretion flow model, which produces the continuum 
spectrum. On the other hand, in terms of the observed QPO phenomena, the origin of the QPO is still under debate, although 
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some theoretical models of the QPO like the disk oscillation model have been proposed. In terms of X-ray polarization of the 
accretion disk, the observed polarization features highly depend on the corona models. However, we do not fully understand the 
production mechanism and physical properties like the spatial distribution of the corona. Therefore, to understand the observational 
signatures containing the information of the spacetime, the physics of the accretion flows and outflows (and/or jets) around the black 
holes should be deeply investigated and understood especially by means of magnetohydrodynamic simulations with effects of the 
radiation fields. Actually, recently magnetohydrodynamic simulations with effects of the radiation are attempted (e.g. De Villiers 
I2OO81 Farris, Liu & Shapiro 2008 Moscibrodzka et al. i2009j . In the near future, the radiation fields around the black hole spacetime 
will be exactly solved by the dynamical simulations with eff'ects of the anisotropy of the radiation field considered in this paper. 
Because the equation given in this study provide the exact results of the contours of the black hole shadow and the photon trapped 
ratio, the equation given in this study can be used for the code tests in the developments of simulation code of the radiation field. 

We have the appearance of the black hole (the black hole shadow) seen by the observer located near the black hole. In summary 
the calculation method for the photon trapped ratio presented in this study can be extended to the astropysically more plausible 
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Fig. 4. Lapse function a representing effects of the gravitational redshift for the black holes with spins of a» = {top) and 0.999 
{bottom). We draw the contours of the lapse function for a =0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 {solid lines) and the event 
horizon {dotted lines). We can see the contours for a - 0.9 around the regions near {x/M,z/M) - (0.75, 0.75) in both panels. At 
the event horizon the value of a becomes zero. 



cases. Although in this study we assume the LNRF, the calculation in this study can be extended to other reference frames moving 
with astrophysically frequently considered velocity fields like corotating or counterrotating Keplerian disks. Actually, in the past 
studies, other frames like a static frame, a circular geodesic frame (or a Keplerian frame) and a radially falling frame considered 
(e.g. Takahashi lIOOTat Schee & Stuchlfk !2009al |2009bl l in addition to the LNRF. Especially Schee & Stuchlfk (2009a) gave the 
equation and contours of the black hole silhouettes by using these frames. In addition to the contours, the equation of the photon 
trapping ratio presented in this study [Eqs. ( |36] | and ( |37] |1 can be straightforwardly extended to a more general velocity field. 
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Appendix A: Derivation of the equation (|57t 

First, we give the analytic expression for the partial derivative d4>ldrs appeared in Eq. ( l36l l. Both rj and f are the function of 
calculated as 

~^ {r]-3r] + alr, + al), ^ (,3 _ 5,2 ^ 9,^ _ 4^2^^ (A.l) 

a,(rs - 1) aiirs - lY 

Based on this, the derivatives of these with respect to r, is analytically calculated as 

idrX) = —r^^ir] - 3r? + 3r, - aj), {d.^i^) = "t?^""' ~ " 3'-? + 3r, - al). (A.2) 
From Eq. ( IA.2l i. we know the relation of 

{dr^rf) = f). (A.3) 
a^Vs - 1) 

On the other hand, the partial derivative of 0, with respect to is given as 

(5, ©,) = -^{d,X) + (dr,n), (d„ V©!) = (A.4) 
tan^e 2V07 

If we introduce 

piB) _ _ „(<!>) _ _ p(r) 

n*-^'* s ——= sine COS n*-'' = — — = sin0sin(/>, n*"* s — =cos0, (A.5) 

p(t) p(t) p(i) 

and 

„('^).^^=tan^=-^(^), (A.6) 
nW ^ Ai/2sin6»\ V©/ 
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we have 

d,-^^ = cos^ ^ (5r,n*"'') 



1 



1 +(n»')^ \Ai/2sin0/ ' \V0/ 
= — T 2(77 + al cos^ 6I)(5, - ^(^r ??) ■ 

On the other hand, the factor 1 - cos Q in the integral given in Eq. ( l36b is calculated as 



1 - COS0 = 1 + 



By using the analytic expressions for c?,^0 and 1 - cos 6, the solid angle S is calculated as 

A'^Vsine r^,^ , 2___2^^,^ . / ■5/- \ -v^ 



S = 2j^' t/r 

= 20"^ sin I drs—= 

Jrr VaCAsin^e©. +pV2) 



2 V0r(A sin^ 6@, + p4^2) 



'"?('".? - 3) ^ 2 2 , 

^ - rj- cos ( 

fl.(''.s - 1) 



-4p^ sin e dr,, 



c 



(ri - 3ri + 3r, - af)[r,^(r, - 3)^ - a,{r, - l)(/7 + ai cos' 0)] (^^^^ ^ s. 



lir, - 1)3(A sin^ 610. + p4^2) ^ 



1 



where the analytic expressions for (, rj, 0, and /?, are given above. 



